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' Abstract 

We examine the string cosmology equations with a dilaton potential in the context of the Pre- 
CN \ Big Bang Scenario with the desired scale factor duality, and give a generic algorithm for obtaining 

solutions with appropriate evolutionary properties. This enables us to find pre-big bang type solutions 
with suitable dilaton behaviour that are regular at t — 0, thereby solving the graceful exit problem. 
However to avoid fine tuning of initial data, an 'exotic' equation of state is needed that relates the 
— ■ ■ fluid properties to the dilaton field. We discuss why such an equation of state should be required for 

f ^ ' reliable dilaton behaviour at late times. 

<N ; 

q\ 1 1 Introduction 

On ■ 

In this paper, we investigate the equations of string cosmology (|] in the string frame, allowing for a 
q-( dilaton potential V(<f>). The Pre-Big Bang Scenario is motivated by the search for cosmological solutions 

with an a(t) — > l/a(— t) symmetry in the scale factor a(t), which implements an analogue of the T-duality 
symmetry of M-theory. However one must distinguish between symmetries of the equations and those 
of their solutions. We look at cases in which the equations have such a scale factor symmetry, when 
solutions may or may not exhibit the same symmetry, and at cases in which the solutions obey the scale 
factor symmetry, even if the equations do not. In the latter case we obtain some solutions that seem 
to have most of the properties desired in the Pre-Big Bang scenario, in that they have the desired scale 
factor symmetry, the desired evolution of the dilaton field, and continuity at t — of a(t), (j>(t), <j>(t), and 
the Hubble parameter H(t) = a(t)/a(t) (but allowing a discontinuity in H(t) and <p(t) there, implying 
a corresponding discontinuity in dV/d(f>), thus providing a solution to the graceful exit problem & El. 
However, to obtain the desired dilaton behaviour at recent times, we need to employ an 'exotic' equation 
of state as discussed below. 

There are 'no-go theorems' that exclude such regular transitions in the presence of a perfect fluid 
and Kalb-Ramond sources. A 'lowest order' Einstein frame analysis by (^] discusses graceful exit in 
generalized phase-space, and derives a set of necessary conditions for transition from a classical dilaton- 
driven inflationary pre-big bang phase to a radiation-dominated era, joined at t = in a Planck epoch 
of maximal finite curvature Hit). They show that a successful exit requires violation of the null energy 
condition (NEC). Classical sources tend to obey NEC, but various new kinds of effective sources generating 
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non-singular evolution have been considered that do not. Thus failing invocation of higher order curvature 
terms, some kind of exotic behaviour of matter is necessary in order to obtain a graceful exit from the 
pre-big bang phase. 

In this paper we follow Gasperini et al |j| by working in the string frame. The relation to the Einstein 
frame is left for another paper. It should be made clear from the start that our solutions are rather 
special in the spectrum of pre-big bang models; those we concentrate on in the main show an exact scale 
factor duality in the solutions, and thus we do not consider here the more exciting possibility of a phase 
of early kinetic-dilaton dominated inflation which leads to an early phase which is not radiation dual 
but is genuinely stringy inflationary vacuum. Nevertheless the set of solutions investigated here help to 
understand the spectrum of possibilities available within the broad Pre-Big Bang set of ideas. 



2 String Cosmology Equations 

One can determine the general equations of string cosmology by extremizing the lowest order effective 
action of dilaton gravity, 



R+(V^f-^H 2 + V{cp) 



J d d+1 x^\g~\L m , (1) 



where </> is the scalar dilaton, H = dB (antisymmetric tensor field strength), V(</>) is the dilaton potential, 
A s is the fundamental string length scale, and L m is the Langrangian density of other matter sources. 
To derive string cosmology equations for the d = 3, homogeneous, isotropic, conformally flat background 
we will follow Gasperini Q| in assuming B = 0, a perfect fluid minimally coupled to the dilaton, and 
a Bianchi I type metric (see Appendix C of Jl| for details). Unlike Gasperini we assume V(<fi) ^ 0, to 
obtain the string cosmology equations in the following canonical form: 

+ (2) 

4>=-3Hi>+^-V-V' + eH^-^, (4) 

where V = ^jr. When combined, these imply the standard energy conservation equation: 

p=-3H(p + p). (5) 

In a relationship analogous to that between the classical Friedmann equation and Raychaudhuri equation, 

Eq.(2) is the first integral of eq.^ provided that eq.^j and eq.^j hold. (6) 

These four equations will be the basis for the analysis in this paper. 

One of the primary motivations for the pre-big bang scenario || is that when V(<p) — 0, these equations 
are invariant under the following transformation: 

a{t) -> d{t) = a-\t) (7) 

provided that the dilaton transforms as <^> — » (f> = <f> — 6 In a and the energy density and pressure as 
p — * p' = a 6 p, p — > p' = —pa~ 6 . Thus if a(t) is a solution, so is d(t) for suitable <fi, p,p. Since the string 
cosmology equations are also invariant under time reversal symmetry, 
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a(t) -> a(t) = a(-t) 



(8) 



the deceleration associated with standard post-big bang cosmology can be associated with an accelerated 
evolution prior to the big bang by the generalized transformation 

a(t) -> a(t) = a -1 (-t). (9) 

where a(t) is a solution for suitable 4>,p,p because a(t) is. The solution has T-duality symmetry if for 
each t, 



a{t)=a(t)=a- 1 (-t). (10) 

However, if we assume V(4>) ^ as in eqs.((^-^|), then in general the equations are not invariant under 
the symmetry eq.( [h]) even if the solutions are. We will look at both cases in what follows, but generically 
allowing a potential that does not preserve the symmetry. Note that if we assume matter with the same 
equation of state before and after t = 0, then the matter equations also will not be invariant under the 
scale factor symmetry. One has to decide what is more physically meaningful: matter with a universal 
equation of state applicable at all times, or that has a discontinuous equation of state that preserves this 
symmetry. In what follows, we adopt the first option. We return to discuss this choice in the conclusion. 



2.1 Flat Dilaton Potential with Exotic Equation of State 

To obtain equations of motion preserving the scale factor symmetry eq. (|Io|) , we assume the simplest 
potential, namely a flat potential: 

V{<j>) = K (11) 

where k is a constant, and then investigate the behaviour of the universe. In order to reliably obtain 
proper limiting behaviour of the dilaton, we assume that the equation of state 

P= 3 + l e ~* K (12) 

holds at all times (this choice, which is not invariant under the duality symmetry, is discussed further 
in the following sections). One can immediately see that at late times if <fi — > constant, as we will show 
follows from this choice, then this equation of state simply reduces to radiation plus a constant. 

We are interested in getting satisfactory dynamics for H(t) and <f)(t) , or equivalently for x(t) = 4>- 
To see when this occurs, we manipulate the string cosmology equations (|[^) subject to eqs.( |Tl| , p^ ) to 
obtain the two-dimensional phase space with coordinates (x,i?) governed by the following equations: 

H = ^-2H 2 + ^, (13) 

* = X(X-3H), (14) 

the latter following because of choice (p"2|). Having chosen the constant k, we can set initial conditions 
(xo,Hq) at t = 0, and then extend the solution to positive and negative values of t by use of these 
equations. For k < 0, there are no fixed points in the phase plane, and on every trajectory both H and 
X diverge as \t\ — > oo. For k = 0, i.e. no dilaton potential, there is one fixed point at the origin, but for 
any initial condition (set at t = ), x an< i H will diverge either as you run time forwards or run time 
backwards. 

The interesting dynamics is obtained when n > 0. There are then fixed points at A + : (0, \J k/12) 
(a source), (0, ^/k/12) (a sink), By. (\/3k, \J k/3) (a saddle point), and -B_: (— \/3/c, — \fnjS) (a 
saddle point). In the phase plane depicted in Figure 1 above, we claim the initial conditions in the region 



3 



-3 

-3 



-2-10123 

Figure 1: Phase portrait representing the solution space of equations (HO) with k > 0. 



I bounded by A + , B + and A— and the separatrixes joining them give satisfactory dynamics of both H 
and x which include x ~~ * as \A — > °°i X > f° r au times so (f>(t) is monotonic, H remains finite, and 
a "bounce" occurs that avoids the initial singularity. Since region I is bounded by fixed points that have 
coordinates proportional to a/k, increasing k will give one a larger region of initial conditions that lead to 
a nonsingular universe with proper dilaton dynamics. We can obtain a solution on the boundary of region 
I (evolving along the line joining A_ to A+, which does not lie in I) that is invariant under symmetry (g) 
by setting xo = 0, Ha = at t — 0, but this solution, given explicitly by a(t) = cosi^^^n/St), is not 
invariant under the symmetry (^). A drawback of all these models is that inflation will not stop at t > 0, 
but as discussed in the conclusion, the string cosmology equations derived in section || do not apply to the 
present cosmological regime without modification, so it is possible that a radiation dominated evolution 
started after the time when these equations no longer apply. In any case this gives a specific family of 
solutions where the equations display the desired symmetry ( |l0| ) but the solutions do not - which is not 
very surprising, given the prevalence of broken symmetries in physics. 



3 Obtaining Desired Dynamics From a Dilaton Potential 

In this section, we generalize the method introduced by Ellis and Madsen Q through which they obtain 
a classical scalar potential associated with a specified a(t) in the standard gravitational equations. No 
field has been observed that coincides with a dilaton potential V{4>) , so we assume that it is a freely 
disposable function. We show that by suitable choice of V(4>) one can obtain almost any behaviour for 
a(i), or alternatively for <p(t). We first present an algorithm for determining V((f>) from a desired a(t) 
or a desired <j>(t), and then present an analytically smooth solvable example. This solution illustrates 
our main point, but has little physical relevance (although it does satisfy the symmetry (0)). In the 
following section we use these methods to obtain two solutions that resemble the standard "pre-big bang 
scenario", but with continuity of a(t) and H(t) and with satisfactory dynamics of <fi(t). The associated 
dilaton potentials are ad hoc because they are derived from the desired behaviour of the universe, rather 
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than from a field theory model; as discussed in many inflationary and quintessence studies, see e.g. 



3.1 The Algorithms 

We proceed by providing the following general algorithm for determining a dilaton potential V(<p) that 
produces a desired a(t): 

1) Specify a desired monotonic function for the scale parameter a(t) , consequently determining H(t) 
and H(t), 

2) Choose an equation of state and solve for p{a) from eq.dqn; as a(t) is known, this determines p(t). 

3) Eliminate V and V from eq.(|J) by use of eqns.(|2|) and (||) to obtain a differential equation relating 
Hit), 4>{t), p(t), and their time derivatives. 

4) Solve the equation obtained in step 3) for (f>(t). 

5) Substitute the now known functions 4>(t), p(t), and H(t) and a(t) into the rearranged version of 

eq-(|) 

V(t) = 6H 2 -e' l 'p-6H<j) + <j) 2 (15) 

to obtain V(t). 

6) Invert (f>(t) to obtain t(<p), and 

7) Transform V(t) as follows: V(t) = V(t ((/))) => V(4>) . This is possible for each range of t on which 
<f>(t) is monotonic (if it is not monotonic on some range of t, in general V(<j)) will not be well-defined 
because it will not be single valued for the corresponding values of (f) . 

Thus, provided <p(t) determined from step 3) is monotonic, we find a V((f>) that corresponds to a given 
monotonic function ait). Because we have now satisfied eqs.(|^), (||) and the equation obtained in step 
3), the latter depending essentially on eq.( it follows from statement that eq.(pj) will be true also, 
so we have satisfied all the equations of the theory for this matter description (c.f. ||); hence we have a 
solution of the desired form. 

Alternatively, we can give an algorithm for determining a dilaton potential V(<f)) that produces a 
desired dilaton evolution <j)(t )Q by proceeding in the same way as above, except for minor changes: replace 
step 1) by 

1') specify the desired monotonic function for the dilaton, <j>it), 
in step 2), leave p in the form p(a), and replace step 4) by 

4') solve the equation obtained in step 3) for ait) (or for Hit)). 
The rest of the algorithm is as before. 

Finally, note that we can carry out these procedures piecewise: for example we can specify ait) for 
some range of t and <p(t) for some adjoining range of t, or different behaviours for a(i) for different ranges 
of t , then join the solutions together, ensuring that ait), Hit), 4>{t) and x(i) are continuous where these 
ranges meet. 



3.2 Exponential Scale Factor Behaviour with No Matter 

To demonstrate the procedure, we give a simple analytically solvable example with a pure scalar field, 
i.e. p = p = 0. Consider an exponential expansion as in classical inflation, 

a(t) = e wf H =io, H = 0, (16) 

where w is a positive constant. This solution has the desired symmetry ( [l0| ). 
In this case the differential equation for (f>{t) takes the form: 

4> = Hd> (17) 

1 If the equation of state is a function of V or V ', then you will have to eliminate these quantities using eqs.(^) and 
) before solving eq. (g). 

2 It is important to note that one has freedom to choose only a(t) or <f>(t), not both. 



■5 



Using eq.Jrij) and eq.([l5|) we obtain 



^(i) = ^ + -(e t0 *-l), 
w 



a monotonic function as required, and 

V{t) = 6w 2 - 6<j> we wt + • 

After inserting the inverted eq.(|l8|), 



into eq.(19), one obtains 



t{<f) = - log 
w 



V(<j>) = w 2 



-y- ( <f> - (po 



? e 2wt . 



<f>o 
w 



3-d> + ^) 2 -3w 2 
w 



(18) 



(19) 



(20) 



(21) 



which is simply a quadratic potential plus a constant. Clearly the behaviour for <f>(t) is unphysical since 
<f>(t) — > oo instead of asymptoting to a constant. However, this gives a transparent example where even 
though the scale factor symmetry (ffoh is broken in the equations (because V{4>) is not constant), the 
solution obeys that symmetry. 



4 "Pre-big Bang" Behaviour 

In this section we try to use the methods just explained to obtain solutions that resemble the "pre-big 
bang scenario" but with satisfactory dynamics of 4>{t) and a continuous transition from the pre-big bang 
to post-big bang phases. In these examples, we seek solutions that evolve from a string perturbative 
vacuum, i.e. H — ► and — ► (no interactions), to the present scenario where e*, which acts as the 
coupling constant, asymptotes to a constant. We will assume the following behaviour of the universe: 

a(t) = (t+l)*, t>0^H(t)=^- (22) 

determining a(t) for t > 0, and by the symmetry jl0| ) 

a(t) = (-t + I)"* , * < =► H(t) = J + (23) 

determining a(t) for t < 0. Both a(t) and H(t) are continuous at t — with a(0) = 1, H(0) = 1/2, but 
H(t) is not continuous there. 

This behaviour, which is essentially radiation dominated evolution of the universe for positive times 
and power-law inflation for negative times, is motivated by the "pre-big bang" scenario introduced in 
||, and exactly obeys the scale factor symmetry (fio]). Note that we have shifted the origin of time 
in each branch from that customarily used, in order to get a smooth evolution through t = 0; this of 
course makes no difference to the desired physical behaviour, for we can choose the origin of time to be 
wherever we want (and the equations are invariant under time translation t — ► t' = t + c). Although the 
power law inflation ends with the scale- factor value a(0) = 1, required by continuity together with the 
symmetry (|Io|), the solution has sufficient inflation for any purpose because it involves an infinite number 
of e-foldings (it starts with the asymptotical value a = as t — > — oo). 
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4.1 Pre-Big Bang behaviour with radiation equation of state 

First we assume the radiation equation of state holds at all times, that is, 

P=f, (24) 

which, using eqs. (§ and (HHI, implies 

p(±t)=p (±t + l)* 2 (25) 

where po is a positive constant and '+£' refers to the post-big bang era, '— V to the pre-big bang era. 
Notice that both p and p are continuous at t = 0. 
The equation for now takes the form 

'^)=\e (t 'p + H(t) + 2H. (26) 

Substituting in eqs. ( p2| ) and (p5[), we could not find an analytical solution to e q.(|26|), so we investigate 
the three dimensional phase space with coordinates (t, <f>, x), given from eqs. ( ]26| , |22| , |23| , p5| ) by 

^r^ ±t+ v T2+ m^ T WTTr> (27) 

where the top sign holds for t > and the bottom sign for t < 0. We can set initial data at t = 0, and 
then investigate the phase plane orbits as we run the trajectory forward and backwards in time in such a 
way that x an d <f> are continuous through t = 0. Then x is discontinuous there, but we have no problem 
in joining the solutions for t > and t < 0. 

For i > 0, there is an exceptional integral curve j(t) given by (t, <f>o,0), where cj)Q = hi(-^); this is 
the only integral curve with a fixed value of <f> and x- Note that setting <pQ and xo determines the initial 
point in the phase space, and specifying po determines the location of this exceptional curve. In the 
2-dimensional sub-spaces t = const with coordinates (4>, x), the curve 7(t) has coordinates (<fio,0) for all 
t, and represents a set of saddle points parametrised by t. To get exactly the desired dilaton dynamics 
in the future (x > , — > constant =>x~^0ast^oo), one must restrict the initial conditions 
(</>o>Xo) to start precisely on the stable branch of these saddle points, which intersects the surface t = 
in a single curve (0, 7+(x), x) passing through the exceptional point 70 = (0, cj>o, 0) (for more details, see 
Appendix A). However there is actually slightly more freedom than this in finding physically relevant 
initial conditions because if a trajectory starts close enough to the stable branch (but not exactly on 
it), then the trajectory will stay close to the fixed point for an arbitrarily long period of time before <f> 
and x diverge, and this may suffice for physical purposes even if the solution eventually diverges (cf. the 
discussion of intermediate isotropisation in Q|). Nevertheless, the physically relevant set of solutions is 
very unstable and requires very precise fine tuning, in order to obtain the desired dilaton dynamics, lying 
in a small open neighbourhood T> + of the curve (j)Q — 7+(xo) i n the initial data set at t = 0. Indeed 
we have found it very difficult to obtain numerical solutions with the desired behaviour because of this 
instability. 

For t < there are no points with a fixed value of 4> and x (because we assume po > 0). To get the 
desired dilaton dynamics in the past (v > 0, — * as t — > —00) one must further restrict the initial 
conditions, the problem being that eq.(|27| ) is an inhomogeneous equation for x with a time- varying source 
function (albeit a source function that decays away as t — > ±00). We can obtain the desired behaviour 
if 2/0 = fe^Vo "C 1, i.e. (f>o <C \n(^) (details are given in Appendix A). This is a sufficient condition; 
there will be a wider domain X>_ of initial data at t = 0, containing this set, that will ensure that at early 
enough times the desired behaviour is attained. 

To get a satisfactory solution for all time, for a given choice of po, one must set the initial conditions to 
lie in both T> + and X>_, so the crucial issue is whether they intersect or not. We have not attained finality 
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on this point. It may be that the 'no-go' theorems with a potential || imply they do not intersect, but 
this implication is not entirely clear, as the conditions of those theorems may not correspond precisely 
to the conditions we contemplate here. If they do intersect, we can attain the desired behaviour x —* 
and — > when time runs backwards as well as — > const as time runs forward and in principle, 
one can obtain a continuous V(4>) associated with the unstable solution described above because every 
function is continuous on the right hand side of eq. (|l5| ) . Furthermore, <p(t) is monotonic and continuous, 
and therefore invertible, so one can complete Step 6 of the algorithm set out in section 3.1. However 
attaining such solutions will require extreme fine-tuning of the initial data, and this is very difficult to 
do because one does not know where the stable branch of the saddle point intersects t = 0. Thus if such 
solutions do exist, the extreme fine tuning required for their initial data make them seem impracticable 
as cosmologies despite their other desirable properties. 



4.2 "Pre-big Bang" Behaviour with Exotic Equation of State 

Finally, we assume the identical "pre-big bang" behaviour of the last example (eqs.(22 2^)), but we obtain 
a stable solution with a different equation of state. The instability in the last example arises because of 
our choice of the equation of state, as can be seen by inspection of eq. (0) , which we write now as 



where 



X = -3H X + X 2 



(28) 



(29) 



As mentioned before, we want to obtain — > constant, i.e. x ~~ * 0, at late times, which implies 
[3 — > in eq. (|28|) . If we choose the radiation equation of state as in the last example, (eq.(24)), then 
= —V — V . Therefore, requiring (3 — > as t ~ > oo puts a heavy restriction on the dilaton potential, 
namely V — > e^^ at late times. Consequently, there is a fine-tuning problem if you use the radiation 
equation of state. 

In the present example, we assume (3 = for all times, which from eq.( |2^) demands the exotic 
equation of state, 



P=l + l e-<»{V+V'), 



at all times (note that eq.(|l^) is the special case resulting when V — 0). Using this equation of state 
implies, 



Pit) 



2He-' t '(12Hi> + 6H - 3<p 2 )dt 



(31) 



which allows the density to go through zero and become negative. We discuss this equation further in 
the Conclusion. 

The differential equation that relates H{t) to <j>(t) is simply eq.( |||) with (3 = 0, 



<f> = -3H(f> + (l) 2 , 

For arbitrary a(t), this can be solved (with a = 1 and xo = <A)) by 

exp(0 o - 0(i)) = 1 - xo / a- 3 {t)dt. 

Jo 

For the specific case given by eqs.(p2|,p3|) we obtain from this the analytical solution 



In 



l-2xo[l-{l + t)---] 



(32) 



(33) 



(34) 



S 



for t > and 



(f>(t) = +(f> - In 



2xo 



[i-(i-m 



for t < 0. Inverting eq.(p4) we obtain 



(35) 



2 2(^-00) 



[(1 - 2 Xo )e+-+° - 1] 



- 1 



(36) 



and inverting eq.(p5ft we obtain 



t(4>) = l 



5e« 



5 + 2xo 



2xo 



-, 2/5 



(37) 



Now we can solve eq.(|3l|) to obtain p(a) and so p(t) (see Appendix B for one particular case), and sub- 
stitute our results into eq. (|l5|) to obtain the dilaton potential V(4>) that is associated with our specified 
"pre-big bang" behaviour. This is straightforward but tedious, and results in very complex analytic 
expressions (the real complexity coming through the expressions for p(t) that occur consequent on the 
choice of the exotic equation of state). Rather than giving these analytic expressions, we give a graph of 
the potential for one particular case below. 

To discuss the relevant initial conditions, it is instructive to look at the phase plane (Figure 2 above) 
with coordinates (t, \ \ where x = <t> is governed by the equation: 



X : 



2(±t+ 1) 



X + X 



(38) 



where we again use + to represent t > and — to represent t < 0. One can easily see that \ — 
(=>■ x = 0) is an attractor, and represents a physically uninteresting solution with <j) — const. Also 
X = 2 (±t+i) * s a nuhcline, characterising the other points where x = 0. This curve starts at (0, |) and 
drops symmetrically away to zero as t — > ±oo. Now we can solve eq. (|38|) analytically for t > 0, finding 



1 



X 



where C± 



2(t+i)(l + C+V*+T) 

1) is positive iff xo < 1/2- 



(39) 



^2^- — 1) is positive iff xo < 1/2- The separatrix between the solutions that diverge and 
those that go asymptotically to zero as t — > oo is the special solution with C+ = which goes through 
(0, i), that is, 



1 



X 



(40) 



2(t + l) 

which itself goes to zero as t — > oo. If we specify the initial conditions at t = such that (j>o is free and 
< Xo < \ C+ > 0, then as we run the trajectories forward in time x ~* 0. In this case, for large 
positive values of t, eq.(|39|) will be approximately 



1 



X = 



2C+ tV2i 



> 



(41) 



(note that <j)(t) is monotonic for t > because x > on these trajectories). Let T + be such that eq.(fi"l|) 
is valid for all t > T + > 0. Then for t > X+, 



t + 2C + t 



■dt + <j)r + — 



C+V2 



-1/2 



- t 



-l/2i 



(42) 
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-10 -6 -2 2 6 10 

Figure 2: Phase portrate representing the solution space of equation 



Thus as t — > oo, for all x : 0(0 ~~ * a constant value, say <j>oo, and exp0(i) — > exp(0oo). (Note that it is 
essential to check this result even though X — > 0, cf. the discussion below of what happens as t — > — oo). 
If we specify the initial conditions at t = such that 0o is free and | < xo ^ C+ < 0, as we run the 

(2xo)!^i i n 



trajectories forward in time then x ~ *• °° as t — > to given by 1 + C+\/to + 1 
this case for large values of X , eq.(Bq) can be approximated as follows: 



0, that is to 



(2X0-1) 2 ' 



(43) 



2(t + l) 

The solution diverges as t — > to and the approximation eq.(^) never applies. This behaviour conforms 
to that implied by eq.(^3), and may be seen clearly on the phase plane. 

If we run the trajectories backward in time, starting from initial data with xo > 0, they will cross the 
nullcline and then drop to zero, never becoming negative because x = is an exceptional solution of the 
equations. Then <j>(t) is monotonic for t < also because x > on these trajectories. Solving eq. 
analytically for t < gives 



X 



(t - %W-t + 1 



(44) 



where C_ 



2(t 2 -2t+lW-t + l + C- 

'2! 1). This expression goes to zero for all C_ > — 1, corresponding to x > (note 

that it does not matter if C_ is positive or negative). For large negative t its value, for all C_, will be 
approximately 

x = d0/dt~-A. (45) 
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3.5 




1 1 1 1 1 1 1 1 1 1 1 

-10 -8 -6 -4 -2 2 4 6 8 10 

Figure 3: The evolution of the scalefactor a(t) as a function of time t, with a(0) = 1, over the time 
interval [—10, 10]. For negative times t < 0, there is power-law inflation, i > 0, followed by a radiation 
dominated phase of expansion for positive time t > 0. 

Let T_ be such that eq.© is valid for t < T_ < 0. Then for t < T_, 

0(t) ~ | (-^)* + ^_ = | In(^) + 0t_ - ex P 0(i) cx {^-f 2 (46) 

Thus as i — + — oo, for all % , — > — oo even though x ~~ * 0, and exp0(i) — ► 0, which is the dynamics 
we desire and indeed is indicated already by eq.({35|). The value of C_ corresponding to the separatrix 
cq.(^) is C_ = 4, which does not give any special behaviour for t < 0. 

Typical results of the integrations for this case are given in Figures 3-5 below. 

In summary, one gets a stable solution for < xo < 1/2, as one can see from the phase plane, with 
good "pre-big bang" behaviour and the desired dynamics for <p(t) for both large and small t. The shape 
of the potential is a bit unusual, but results directly from the specific requested 'pre-big bang' evolution 
eqs.( |2^p^ ) and the chosen initial conditions. Smoothing out that behaviour at t = 0, so that the solution 
departs from the 'radiation' form eq.(|2^) at very early times while preserving the symmetry (|l0|), will 
result in a smoothed out potential V(<j)) ; we can choose a(t) in this way so that H (t) and hence V{(f>) are 
continuous at t = 0. Initial conditions can be set so that the matter has the desired late time behaviour: 
p/ p — > 1/3, p — * 0; however it then has unusual behaviour at early times in that both p and h = p + p go 
negative for some values of t < 0. It is unclear if this should be regarded as a serious defect of the model or 
not, remembering that with the unusual equation of state adopted, the properties of matter are different 
than usual, and in particular the speed of sound will no longer be given by the usual expression. This 
needs further investigation. What is clear is that these solutions are not physically reliable as t — > +oo 
(see below), and they will have to be joined on to some other solution to give an adequate model of the 
universe with ordinary matter behaviour at late times. However, as discussed below, that problem occurs 
in the entire family of pre-big bang models, and so is not restricted to the models considered here. 
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Figure 4: The function exp((/>(t)) as a function of time t, with a(0) = 1, (f)(0) 
exp (4>(t)) increases monotonically from at time t = — oo to 2 at t = +oo. 



and x(0) = 0.25. 




Figure 5: The dilaton potential V((f>) as a function of time <j>. We assume that a(0) = 1, (f)(0) = and 
x(0) = 0.25 and take the density p(t) to have value p(oo) = at time t = oo. The potential V((f>) is 
continuous at all times, but non-differentiable at <j> = 0. For <fi — * — oo, 1^(<^) is asymptotically zero. To 
the right of = 0, the potential starts at V w —0.005 and goes to zero from below as phi goes to In 2, 
then increases to +oo as <fr — ► oo. Around = In 2, both V(<p) and its gradient V'(<j>) are zero. As time 
t —> +oo, the dilaton field asymptotes to a constant value of In 2 in our model. The dilaton potential 
V(<f>) approximates a fixed value of as <ft — ► In 2 asymptotically for large positive times. 
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5 Discussion 



We have given examples making very clear the distinction between the equations and the solution having 
the desired 'pre-big bang' symmetry. We have given a broad method of attaining desired string cosmology 
solutions when there is a dilaton potential V not equal to zero, and used it to obtain 'Pre-Big Bang' 
solutions that seem to have close to the desired properties. In the first case considered, choice of the exact 
radiation equation of state at all times leads to a very unstable situation where extreme fine-tuning 
of initial conditions is required to attain the desired results, and indeed there may be no initial data 
leading to the desired behaviour in both the forward and backwards directions of time. In the second 
case we impose an 'exotic' equation of state (|o|) that links the fluid behaviour to the potential in a way 
that generalises the perfect fluid equation of state, and we obtain solutions of the desired type without 
the need for fine tuning the initial data set at t = 0. 

This equation of state looks strange, and the resulting matter behaviour is certainly unusual, but we 
have no solid handle to use in restricting equations of state in this early era; and we suggest that it is 
essential to choose such an equation if one wants the solution to reliably tend to the 'classical' form at 
late times. This is because of the form of the equation for (j>; if we do not set (3 = 0, where (3 is defined 
by eq.(p9|) then almost always that desired classical state will not be attained, because of eq. (p8|) ; but 
setting (3 = 0, which leads to the desired behaviour, leads immediately to our 'exotic' equation of state. 
Insofar as that equation of state and resulting behaviour is unsatisfactory, this indicates that there is a 
problem with the form of the equation for (f>, which comes directly from the standard variational principle 
employed in the context of the pre-big bang scenario. The remedy probably lies in finding other scenarios 
with alternative forms of the variational principle, leading to other equations for </>. 

This is also indicated because the present form of the equations does not accommodate ordinary 
matter, the point being that the above analysis applies even if there is no dilaton potential. Suppose 
V = 0; then eq.(|2| ) remains true, but now 

/? = e*(f-f), (47) 

so a reliable approach of the dilaton to a classical solution at late times, requiring [3 = 0, demands the 
radiation equation of state ( pi| ) ; a baryon dominated epoch is not allowed^. This is usually dealt with by 
stating that eqs.(||-||) don't apply at late times in the history of the universe - a different set of equations 
are to be used then, and the solutions for early times obtained from cqs.( must be suitably joined 
on to that late time evolution. However given the vision of M-theory as representing the fundamental 
theory of gravity, it should be able to describe that epoch too; this apparently requires some modified 
scenario and associated variational principle (note that although we have discussed the issue in the string 
frame, it also arises in essentially the same form in the Einstein frame). In any case, whether one accepts 
this argument or not, given the standard variational principle and equations, we argue that the 'exotic' 
equation of state implied by setting f3 = is necessary to give the desired behaviour; when adopted, it 
enables obtaining that behaviour reliably (i.e it eliminates the need for extreme fine-tuning of data set 
at t = 0). 

However one should note here that we have perhaps been somewhat extreme in imposing this equation 
of state at all times. It is only really needed, on our approach, near the time of the turnaround, and one 
could obtain far more general behaviours by modifying what we have here in that light; what is required 
is that the quantity (3 must go to zero in the period when the dilaton is stabilised. It has also been pointed 
out to us that it is not clear why the deviation from its vanishing point should be absorbed completely in 
the pressure, and then promoted into the conservation equation; other models of the transition [ll], [l2| 
successfully stabilise the dilaton at late times without this requirement, with suggestions for classical and 

3 Although of course by the algorithm given above, we can simulate a matter dominated phase by suitable choice of the 
potential V. 
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quantum corrections in the effective action taking the place of the exotic fluid. Hence our proposal must 
just be seen as one of a range of possibilities in this regard. 

Because we have not made the usual separation of our solution into a '+' and a '-' branch, it is 
not immediately clear why these solutions are not ruled out by the 'no-go' theorems involving a dilaton 
potential || ; this is presumably because those theorems exclude fluids with the equation of state we have 
assumed. We also have not examined the relation of these string-frame solutions to the corresponding 
Einstein-frame versions. These issues await investigation. 

We thank M Gasperini, A Coley, R Tavakol, and the referee for helpful comments, and particularly 
J Lidsey for helpful discussions. DR wishes to thank Elaine Kuok for her patience in checking many of 
the calculations in this paper. We thank the NRF (South Africa) and Queen Mary College, London, for 
financial support. 
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Appendix A: Pre-Big Bang Evolution for Radiation 

For given po, it is convenient to define y = |e^po arL d change variables to (t,y,x)- The equations (| 
for t > become 
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y = ^y^ = frh + WTiy (48) 

In the 2-dimensional sub-spaces t = const with coordinates (y, x), the curve "f{t) has coordinates (1, 0) for 
all t, and represents a set of saddle points parametrised by t. To get exactly the desired dilaton dynamics 
in the future (x > 0, — + constant \ ~~ * as i — > oo), one must restrict the initial conditions (j/0)Xo) 
to start precisely on the stable branch of these saddle points, which intersects the surface t — in a curve 
(0,7+(x),x) passing through the exceptional point 70 = (0, 1,0). One can obtain approximate solutions 
by rewriting the second of eqs.(^) in the form 

x X y- 1 



{i + ty/ 2 ) (t + 1) 5 / 2 

Suppose y is almost constant for t > T+, implying x is close to zero then. Then we can integrate to get 

t> t+ => x = -|yT7 + c+ ' /YT ~ t 

where C' + determines the magnitude of x at time T + The first part decays away as desired, but the second 
part grows with time unless C+ = 0; this is the fine-tuning required to attain the desired behaviour of x- 
To investigate t < 0, it is again convenient to define y — ^e^po and change variables to (t, y, x). The 
equations for t < become 

y = xv, x = y(-t + 1) 2 + + (49) 

implying that \ > for all t < 0; hence \ necessarily decreases at all times in the past. The problem is 
that it can become negative, because x = is not an invariant set of the equation. We want a solution 
where x remains positive for all time so that <f> decreases for all time; this means we need x to go to a 
positive value or zero, but not to become negative, and y to go to zero. . As in the previous case one can 
obtain approximate solutions by rewriting the second of eqs.(E9h in the form 



(x(i-t)^)- = H+ 1 i)3/2 (i + ,H + i) 4 ). 

Suppose 



y(-t + l) 4 < 1 for t < T_. (50) 
Then we can ignore the second term on the right and integrate to get 



t<T^ x = — + —, y = y ^- w eM-2C-VT^t) 

where C_, yo represent the magnitude of x, y at time T_ . This decays away as desired, and consistently 
preserves the inequality (|5^) for all earlier times because the exponential always dominates the power 
law terms. The question then is whether for suitable initial conditions we can attain this inequality at 
some time T_, requiring y(T_) -C (1 - T_)~ 4 . We can satisfy this with T_ = if y = |e^°p < 1, i- c - 
O « ln(^). 

Appendix B: Density evolution with exotic equation of state 

The 'pre-big bang' evolution (B3,E3) implies H and H in terms of a: 
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t>0: H(a) = H(a) = (51) 
t < : H(a) = 4, = if- (52) 



Assuming the exotic equation of state (|30|) implied by setting j3 = at all times, from (|33j) we find <p in 
terms of a: 

i>0: exp(0(a)) = exp(0 o ) n(1 _ 2x ° )+2xo , (53) 

t<0: exp(0(a)) =exp((/. ) a5/2(1 _1 x 2 o)+ . Xo , (54) 

and from ( |39| , [44"| ) we find x m terms of a: 

*>° : X(a) = Ba(2xo+ gL 2xo) n) , (55) 

xW = - 2Xo+( S o)o ^ - (56) 
A particularly simple case occurs when Xo = 2- Then 

t>0: exp(0(a)) =exp(0 o )^ I , (57) 

i<0: CX p(0(a)) =cxp(0 o )g^^. (58) 



and 



* > : X(o) = 2a . (1 1 +a) , (59) 
t<0: X(a) = - 2{1 % a5/2y (60) 

Now p(t) is determined by d3l|); using the above expressions, for t > and = \ this becomes 



dp 3 3 

da ~ a 5 4a 6 (l + a) 



which can be solved to give 



, . _ 3 9 1 3 3 3,. a 
p( a ) = C + 777* + 7777 + 777 " 73 + 77 + 7 ln (; 



20a 5 16a 4 4a 3 8a 2 4a 4 v l + a y ' 

This implies p(t) -^C+^f-fln2 = C + 0.81764... as t -» 0+ and -> C as i -> oo; hence choosing 
C = 0, p(t) — > 0.81764... as t — > 0+ and p(i) — > as t — > oo. Also pip — > 1/3 as t — > oo. The expression 
for in this case follows on putting this into([l5|) and using ( [36]), (|22|), and the various expressions 

above. Similar (more complicated) expressions can be obtained for t < 0. 



16 



